Non-local double-path Casimir phase shifts in atom interferometers 
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We present a quantum open system theory of atom interferometers evolving in the quantized elec- 
tromagnetic field bounded by an ideal conductor. Our treatment reveals an unprecedented feature 
of matter-wave propagation, namely the appearance of a non-local double-path phase coherence. 
Such a non-local phase arises from the coarse-graining over the quantized electromagnetic field and 
internal atomic degrees of freedom, yielding a non-Hamiltonian evolution of the atomic waves mov- 
ing in presence of correlated quantum dipole and field fluctuations. We develop a diagrammatic 
interpretation of this phase, and estimate it for realistic experimental parameters. 



Thanks to the advances in atom cooling achieved in 
the last decades, atom interferometry has become a field 
of great importance for both basic and applied science, 
enabling, in particular, the realization of extremely accu- 
rate inertial sensors jU fjj . With the advent of the coher- 
ent atomic waves guided on chips [3j , the investigation of 
atom-surface interactions has become a frontier for such 
systems. Already, atom interferometers have been used 
to experimentally probe the van der Waals regime [3]. 
This experimental effort calls for a complete theory of 
atom interferometers in the presence of quantum fluctu- 
ations of the electromagnetic (EM) field. 

In this Letter, we layout such theory for a beam of neu- 
tral atoms, compute the phase shifts, and compare with 
the results obtained by standard techniques suitable for 
atoms driven by conservative forces. First, we develop a 
theory of the atomic phase-shifts taking the effect of field 
and atomic dipole fluctuations separately over each inter- 
ferometer arm. This method, however, neglects quantum 
correlations, mediated by the field, between the atomic 
wavepackets evolving along the separate arms. In or- 
der to capture this new effect, we develop a theory of 
atom interferometers based on the influence functional 
method [5]. We show that non-local double-path phase 
shifts arise which, to our knowledge, have not been pre- 
sented before. Such non-local phase is associated to pairs 
of paths rather than to individual ones. It is absent in 
the standard treatment of matter-wave dynamics with 
conservative forces, which shows that the effect of quan- 
tum vacuum and zero-point dipole fluctuations on atomic 
waves cannot be understood as an effective potential. 

Typical atom interferometers [TJ [5] involve the succes- 
sive propagation of atomic samples through atom-optical 
elements such as atomic beam-splitters or atomic mirrors 
(usually implemented with the help of laser or magnetic 
fields) and in the dark (spatial regions with no applied 
fields) . We consider here the phase arising from the prop- 
agation of atomic waves in the dark region along the two 
arms of the interferometer. We assume that the external 
atomic degrees of freedom (d.o.f.) evolve according to 
a quadratic and separable potential V(r), which is a 
realistic situation in many atom optics experiments. The 



initial atomic state is assumed Gaussian, i.e. |e(£ )) = 
Jd 3 r wp(r,r ,po,w )|r) with wp(r, r , p , w ) = 

n, )= ., s , z (Vy^o^)e" (r '"' )o)2/2,i,o2 '' +4po '' ( ^' ,o)/?l - If the 
atomic wave-packet is sufficiently dilute, the evolution is 
given by the ABCD theorem for atomic waves [B]: 

(r|e(i)) = wp{r, r(t), p(t), w^e^^'P"''-*") (1) 

for t > to. The average atomic position r(i) and 
momentum p(t) follow the classical equations of mo- 
tion associated with the Hamiltonian Ha, with the 
initial conditions r(t ) — r and p(to) = Po The 
ABCD approach, exact for non-interacting samples, 
can be extended to include moderate interactions [7 
and atom lasers [5]. The width vector w(i), known 
analytically, is unimportant for the coming dis- 
cussion. The action, 5(ro,po,i — to), is obtained 
from the external and internal (energy E(t)) atomic 
Hamiltonian, integrated along the classical trajectory 

5(r , po, t - t Q ) = Jldt' (ej? - E(t') - V(r(t>))) . 

We begin with a local description of an atom interfer- 
ometer in presence of quantum dipole and EM fluctua- 
tions. As depicted in the two diagrams of Fig. [TJ we con- 
sider two arms sharing the same origin and corresponding 
to atoms flying between the instants t and tf above a 
metallic plate. One of the arms is parallel to the plate, 
and the other has a small velocity component normal 
to the plate. This simple configuration is sufficient for 
the discussion, which can be easily generalized to more 
complex interferometer geometries. The initial external 
atomic quantum state is a coherent superposition of two 
Gaussian wave-packets of common average position and 
different momenta, i.e. \ip e o) = ^ l e i(^o)) + l e 2 (^o)) with 
(r|efc(io)) = wp(r, r , pofc, w ) as defined above. To ac- 
count for the EM field fluctuations, we now consider the 
full Hamiltonian H — Ha + Hp + Haf with the field 
Hamiltonian Hp = ^ k Huj^ ^a^ak + 1/2^ and the dipole 

interaction Hamiltonian Haf = — |d • E(f Q ). Note that 
the dipole interaction involves three quantum operators, 
since the atomic position f a at which the electric field is 
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taken is also an observable. 

We now apply the standard procedure in atom interfer- 
ometry. Precisely, we assign to each interferometer path 
k a phase corresponding to the average dipole interaction 
energy for an external atomic state \e k ) following this 
path, i.e. = i j£ dt(d ■ E(r a ) + E(r Q ) • d) fc (t). To 
derive this phase, we use linear response theory and treat 
the dipole interaction Hamiltonian Haf as a perturba- 
tion [5]. In the Heisenberg picture, the atomic dipole 
operator then becomes the sum of a free-evolving contri- 
bution and of a contribution induced by the coupling to 
the on-atom electric field. The induced dipole operator 
is expressed as the convolution of the atomic polarizabil- 
ity with the free on-atom electric field playing the role 
of a source. By linearity of the coupling, the dipole and 
the on-atom electric field operators play symmetric roles, 
thus a similar result holds for the latter. 

The phase (j)^ is then expressed as a propagation in- 
tegral involving two contributions, namely the response 
of the atomic dipole to the on-atom electric field fluctua- 
tions and the response of the on-atom electric field to the 
dipole fluctuations. The fluctuations of the dipole and 
electric field are captured by Hadamard Green's func- 
tions, whereas the susceptibilities (polarizability for the 
atom) correspond to retarded Green's functions. We as- 
sume that the width of the atomic wave-packet is small 
compared to the relevant EM field wavelengths, which 
allows us to replace the on-atom electric field operator at 
a given time by the electric field evaluated at the corre- 
sponding classical atomic position (t) on the considered 
path k (electric dipole approximation) . The retarded and 
Hadamard Green's functions read 



(2) 



with the operator 6 = d, E, where corresponds to 
the free-evolving Heisenberg operator. For the dipole op- 
erator O = d, the arguments of the Green's functions 
are two instants (x,x') = (t,f). For the electric-field 
operator O = E, these arguments are two four-vectors 
(x,x') — (r,t,r',t'). Last, we use the isotropy of the 
dipole Green's functions to write the standard Casimir 
phase as: 
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dtdt' 



G?(t,t')TrG%(r k (t),r k (t>)) 



+ G£(f,t')TrGf (r k (t),r k (t')) 



(3) 



with the four- vectors r k [t) — (r k (t),t) and Tr denoting 
the trace of the Green's functions G-' with respect to 

E. ij 1 

its tensorial components. We note that this phase already 
contains non-trivial dynamical Casimir effects [10] . be- 
yond the plain integration of the Casimir potential taken 



at the instantaneous position along the classical trajec- 
tory, particularly in the case of a time-dependent atom- 
surface distance as in the trajectory r 2 (i) shown in Fig. 1. 

We now turn to a rigorous computation of the 
phase with closed-time-path (CTP) integrals [5]. Such 
path integrals involve simultaneously forward and back- 
ward histories of the system. The coarse-graining 
over the quantum field induces a mixing between 
these two histories, which will be associated to a 
pair of interferometer paths. This is at the ori- 
gin of a path entanglement inducing the double-path 
phase discussed below. First, we specify the actions 
for the considered quantum d.o.f., namely 5 b [r a ] = 
ft* (ir*"|(i) — V[r (i)]) for the atomic position r a , 
S F [A»] = (e /4) / d 4 xF^F^ with F^=d^A v - d v A„ 
for the EM field, and the action for the dipolar interac- 
tion SafIA^, d, rj = — J d 4 x J^(x) A fl (x) defined with 
the current J" [d, r ] (a:) = - / dtd^t) /cf 5 (4) (x - r a (t) ). 
The four-dimensional integrals are defined as J d x = 
f* f dt J d 3 r. We have introduced the differential oper- 
ator — dii]Q^ — d^rji^ {tj^v is the Minkowski metric 
with mostly plus signature) relating the electric field to 
the vector potential by contraction: 2£,(x) = k^AAx). 
The following discussion is valid for an arbitrary action 
Sa [d] for the internal atomic d.o.f.. We consider the evo- 
lution of the density matrix p(r a , d, A 11 , t) and trace over 
the field and dipole d.o.f.. The result can be expressed as 
a CTP integral over the position involving an influence 
action Si F [r a ,r' a ] 

J CTP 

where v a and r' a refer to forward and backward histories, 
respectively. We have noted for later convenience the 
CTP integral for a generic d.o.f. X: 



CTP 



VX= I I dX dX' 



x', 



„ , *VX [ f vx'p(x Q ,x' ,t ) 

Xn J XL 







(5) 

with p(to) the initial density matrix. At this level, the in- 
fluence action S/p'[r a ,r^] fully accounts for the effects of 
the field and the internal atomic dynamics on the exter- 
nal atomic d.o.f.. The final density matrix p(tf) contains 
four sharp peaks in the region (r Q ,r^) centered around 
the classical positions (rj(fy-), Tj(t/)) for i,j = 1,2. The 
desired phase shift is obtained from the off-diagonal den- 
sity matrix elements (i.e. i ^ j). The main contribution 
to the path integral comes from the paths in the vicinity 
of the two stationary paths. Thus, the phase difference 
due to quantum fluctuations can be evaluated on the clas- 
sical paths by taking the real part of the influence action 
4>if — |Re [Sif[^2, r i]] [H] (whereas the imaginary part 
represents the decoherence due to the plate [H]). To de- 
rive these classical paths, one can neglect the influence 
action corresponding to the Casimir-Polder interaction, 
which amounts to finding the Casimir phase shifts by 
considering undeflected trajectories. 
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Let us detail the procedure to obtain the influence 
action in the spirit of previous derivations of non- 
equilibrium forces mediated by a quantum field [TUl 
113] . First, we define an intermediate influence action 
Sjp [d, d', r a , r' a ] corresponding to the trace over the EM 
field alone: 



aLaI 



CTP 



xe i(S P [A fl ]+S AF [A ll ,d,v a ]-S P [A' ll ]-SAF[A' li ,d',r' a ]) ,qs 

This CTP integral over the vector potential linearly cou- 
pled to an external current is well-known (see for instance 
Ref. [IT]) 



We have introduced the plain and time-ordered corre- 
lations of a free-evolving single dipole component (all 
diagonal components are equivalent by isotropy) , respec- 
tively <?+, = £(di(iR(i')>d and g[ t , = %{Tdi(t)di(1?)) d . 
tff = iRe^Kn]] and <j>ff = ±Re [S?/[r 2) n]] 
are the single and double path influence phases respec- 
tively. To evaluate the single-path phase (j)fj, we note 
that both electric field Green's functions are real and use 
the general dipole relations Re [gf tt i] — |C??(t, i') and 

Im [g[ t ,] = -i (Gf(t, t') + Gf{t', tj\ . One then obtains 

that <fifj = <pft — (f$, namely the SP influence phase 
coincides exactly with the standard local phase of Eq. Q . 



^[d,d',r a ,r' a ] = 



d 4 xd 4 x' 



J' 1 -{x)Gl iiv {x,x)r + {x) 



+ ~^-(x)Gl^(x,x')J 1/ -(x') 



(7) 



Taking standard conventions, we have introduced the 
semi-sum J + = |j[d,r ] + ^J[d',rJJ and difference 
J~ = J[d,r ] — J[d',r' a ] variables. The vector poten- 

R H 

tial's retarded and Hadamard Green's functions G ,' 

are defined as the electric field Green's functions G^ . . in 
Eq. (pi) with Cartesian coordinates replaced by Lorentz 
indices. 

To obtain the desired influence action S7,F(r a , r^k we 
average the EM influence functional given by Eqs. (6|7 1: 
e iS IF [r a y a ] = ( e jiS%ld,d',r a ,r' a ]} d with ^ ^ denoting 
the time-dependent average over the free-evolving dipole 
d.o.f.. Wc take an approximation of small dipolar cou- 



pling, namely e^ s = 1 + j^S + 0(d 3 ) for the actions SVf 

and Sfp. Thus the influence action reads Sip[r a , rJJ ~ 

(Sjp[d, d', r a ,r' a ])d- The vector potential Green's func- 
tions can be treated as constants when averaging over 
the dipole d.o.f.. In order to express the influence ac- 
tion in terms of EM and dipole correlation functions, 
we expand the currents J in terms of J[d, r a ] and 
J[d',r' a ] in Eq. ^ and integrate over the spatial coor- 
dinates. The influence action receives a single-path (SP) 
and a double-path (DP) contribution, i.e. S , / F [r 2 ,r 1 ] = 
5f|[r 2 , ri ] 



S?/[r 2 ,ri] with 



dtdt' 



la 



gF t ,TrG?(r 2 (t),r 2 (f)) 



gl t ,TTG"(r 2 (t),r 2 (t'))+g[*TrG"(r 1 (t),r 1 (t')) 



S^[r 2 ,r 1 }=- 



dtdt' 



g+ t ,TrG%(r 2 (t), ri (t>)) 
-gt;,TrGZ( ri (t),r 2 (t')) 



'g+ TrGf (ra (t) , n (t) ) + g+* TrGf (n (t) , r 2 (f) ) 



(8) 



On the other hand, the presence of a non-local double- 
path phase <^ p contrasts sharply with the local phase 
obtained by a standard atom interferometric method. In 
order to interpret and evaluate this phase, a first step 
is to derive the trace of the electric field Green's func- 
tions [Eq. Q] in the presence of a perfect infinite plane 
conductor placed at z = 0. The result can be written as 
the sum of a free-space and of a scattered contribution, 
namely TrGj?(r, r' 



G & {r,r 



Q^' S (r, r') with 

E 



£ E [r,r ) 



O(t) d 2 /S(t-\R\/c) 



2ne c 2 dr 2 \ |R 
e(r) d 2 ( Sjr-m/c) 
2tt€ dzdz' \ |Ri| 



(9) 



By translational invariance, the free-space contribu- 
tion Ga'°(t, r') depends only on the relative position 
R = r — r'. The scattered contribution involves the 
image of the dipole, and thus we have introduced 
Ri = r — r{ with the image position r[ = (x 1 , y', — z'). 
Both contributions depend on the time difference 
t = t — if and not on the individual times. A similar 
expression holds for the Hadamard Green's function, 
but the corresponding phase contribution oscillates so 
fast that it can be neglected in practice. 

The influence action of Eq. ^ appears as a sum of 
propagation integrals, which can be turned easily into 
Feynman-like diagrams involving the simultaneous prop- 
agation on distinct paths. Using the Green's functions 
of Eq. @, one obtains the DP phase as the difference 
of two free-space diagrams (a)— (b) plus the difference 
between two scattered diagrams (c) — (d) as depicted in 
Fig. [T] Each diagram involves an advanced time t on 
one path and a retarded time t' on the other path (for 
the free-space diagrams) or on the image of the other 
path (for the scattered diagrams). These paths will be 
denoted respectively the advanced (red arrow) and the 
retarded path (blue arrow). The purple arrow stands for 
the free-space (diagrams (a) and (b)) or for the scattered 
(diagrams (c) and (d)) retarded electric field Green func- 
tions. These diagrams show that quantum dipole flue- 
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FIG. 1: (Color online) Diagrammatic expression of the DP 
phase. The delay t — t' between the advanced and re- 
tarded times corresponds to the propagation of light over 
the distances \ri(t') — r m (t)| (for the diagrams (a,b)) and 
\ r u{t')— i"m(t)| (for the diagrams (c,d)) with I, m — 1,2,1 / m. 



tuations propagate, via the quantum electric field, from 
the retarded path (four-position ri(t') for the free-space 
diagrams, image four-position ruit 1 ) for the scattered di- 
agrams (c,d)) to the advanced path four-position r m (t) 
with m I. Such diagrams express the Casimir inter- 
action between two coherent components of a quantum 
state associated with a single atom, travelling on differ- 
ent arms (free-space diagrams), or on one arm and on the 
image of the other arm (scattered diagrams). 

Let us focus on the scattered diagrams (c,d). Since 
the instants t, t' are constrained by the causality prin- 
ciple (t' < t) and by the speed of light, the atom 
moves on the advanced arm during the propagation time 
t — t' ~ |ri (i) — Y\2(t)\/c from the retarded path to the 
advanced path. The resulting distance is different de- 
pending on whether the path 1 or the path 2 is the ad- 
vanced path, because they correspond to different veloc- 
ities. This yields an asymmetry between diagrams (c) 
and (d): the distance between the retarded image po- 
sition vu(t') and the advanced position r m (t) is indeed 
slightly larger in the left diagram (c) than in the right di- 
agram (d) , thereby yielding a stronger cross-talk between 
the quantum dipole fluctuations on distinct paths in the 
latter. The DP phase is essentially a signature of this 
asymmetry. Note that such asymmetry is reminiscent 
of two basic features of the EM interaction, namely the 
causality and the finite speed of the propagation. Such 
asymmetry arises also in the free-space diagrams (a) and 
(b), precisely the quantum dipole cross-talk is always 
weaker on the diagram attributing the advanced times 
to the faster path. In the short-distance limit considered 
below, this free-space contribution is actually of much 
smaller amplitude, and independent of the distance zq. 
It is nevertheless remarkable that a non-local DP phase 
persists even in empty space. This effect will be discussed 
in detail elsewhere. 

We now give an analytical expression of the double- 



path phase for a realistic experimental system, taking the 
two classical paths shown in Fig. [T] expressed as ri(t) — 
r +U|| (t— t )x and r 2 {t) = r +v^(t—t Q )x+vj_(t—t )z. To 
evaluate the dipole correlation functions gj~ t , , we modeled 
the dipole as an harmonic oscillator with a transition 
frequency luq, which yields the correlation function g^ t , = 

|a(0)a;oe _4 " o ^ _t ' expressed as a function of the static 
atomic polarizability a(0). In the short-distance limit 
z < A , the DP phase reads (T = tj — t ) 

DP _ 3tt fa(0)\ ( 1 I 



3tt 

4^ 



47re r 



z 2 Q (z + v±T)< 



(10) 



where Ao = 2ttc/u!q is the transition wavelength. For 
a short path separation v±T <C z , the phase (j^f be- 
comes proportional to the endpoint separation v±T. For 
a long path separation, i.e. zq <C v±T <C c/ujq, the DP 
phase saturates to a maximal value depending only on 
the atomic polarizability, on the transition wavelength 
and on the altitude z . We have investigated this value 
for a sample of 87 Rb atoms with a static polarizability 
a Rb(0)/(47re ) = 4.72 x 10~ 29 m 3 where the main transi- 
tion has a wavelength Ao — 790 nm. We have considered 
a distance to the plate of z — 20 nm similar to the dis- 
tance used in the experiments of Ref. [4.. Such param- 
eters yield a DP phase of <j>?f = 3.5 x 10" 7 rad. This 
value is beyond the sensitivity given by the state of the 
art atom interferometers but still larger than systematic 
phases considered in atom gravimeters [2]. 

To conclude, we have developed a quantum open sys- 
tem theory of atom interferometers, predicting non-local 
double-path phase shifts in the propagation of atomic 
waves. We have shown that the standard atom-optics 
approach catches only the local phase shifts correspond- 
ing to single-path terms obtained with the influence func- 
tional method. Such non-local DP phase is a consequence 
of the coarse-graining over the EM field, which induces 
cross talks between quantum dipole fluctuations on each 
arm. The DP phase constitutes a one-particle quantum 
interference effect, which exists thanks to the coherence 
of the matter waves propagating in the interferometer. 
We have developed a diagrammatic picture of this effect, 
which occurs even in free space, and can be interpreted as 
an asymmetry between diagrams involving simultaneous 
atomic propagation on distinct paths. This asymmetry 
can be seen as a dynamical relativistic correction of the 
quantum dipole cross-talks between the interferometer 
arms. We have shown that the DP phase shift compares 
to systematics considered for accurate atom interferom- 
eters. Our approach can be extended to multiple-path 
atom interferometers by considering pairs of paths. This 
is to our knowledge the first evidence of a non-local phase 
coherence in atom optics. 
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